1 



LINEAR SECOND-ORDER DIFFERENTIAL EQUATIONS 
FOR BAROTROPIC FRW COSMOLOGIES 

H.C. Rosu 1 :}:, O. Cornejo 1 , M. Reyes 2 , D. Jimenez 2 

1 Dept. of Appl. Math., IPICyT, Apdo Postal 3-74 Tangamanga, San Luis Potosf, MEXICO 
2 Institute de Ffsica, Universidad de Guanajuato, Apdo Postal E-143, Leon, MEXICO 

q '. Abstract 

Simple linear second-order differential equations have been written down for FRW cosmolo- 
gies with barotropic fluids by Faraoni. His results have been extended by Rosu, who employed 
^ ! techniques belonging to nonrelativistic supersymmetry to obtain time-dependent adiabatic in- 
dices. Further extensions are presented here using the known connection between the linear 
second-order differential equations and Dirac-like equations in the same supersymmetric con- 
text. These extensions are equivalent to adding an imaginary part to the adiabatic index which 
is proportional to the mass parameter of the Dirac spinor. The natural physical interpretation 
CN ■ of the imaginary part is related to the particular dissipation and instabilities of the barotropic 
FRW hydrodynamics that are introduced by means of this supersymmetric scheme. 
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O . The barotropic FRW cosmologies obey the following set of differential equations: 

7 - = ^(P + 3 P)> (!) 
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1 - Introduction 



a J 3 a 2 

P = (7"1)P, (3) 

where a is the scale factor of the universe, p and p are the energy density and the pressure, 
respectively, of the perfect fluid of which a classical universe is usually assumed to be made 
of, K — 0, ±1 is the curvature index of the flat, closed, open universe, respectively, and y is 
the constant adiabatic index of the cosmological fluid. Recently, Faraoni [ ] proposed the 
"Riccati route" of solving the system of eqs. (l)-(3) and Rosu used Faraoni's approach to 
introduce a supersymmetric class of cosmological fluids possessing time-dependent adiabatic 
indices [2]. It was claimed that these fluids can provide a simple explanation for a currently 
accelerating universe [3]. 

In this work, we review the supersymmetric factorization methods for barotropic FRW 
cosmologies in section 2. Next, in section 3, we present corresponding Dirac-like (first-order) 
coupled differential equations and their associated second-order differential equations and 
discuss them in a formal way. We end up the work with a short conclusion section. 
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2 - Supersymmetric (factorization) methods 

Combining the equations (l)-(3) and using the conformal time variable rj defined by 
dt = a(r/)dr] one gets the equation 

2 

+ ck = 0. (4) 




where c = |y — 1. The case k = is directly integrable [ ] and will be skipped henceforth. 

One can immediately see that by means of the change of function u = |- the following Riccati 
equation is obtained 

u'+cu 2 + kc = . (5) 
Employing now u = - ^ one gets the very simple second order differential equation 

w" K + kc 2 w k = . (6) 

For k = 1 the solution of the latter is Wj = Wj cos(ct] +d), where d is an arbitrary phase, 

implying a l (rj) = A 1 [cos(cr/ +d)] 1//c , whereas for K = — 1 one gets w_j = W_ 1 sinh(c7]) 

and therefore a_j(r]) = A 1 [sinh(c7])] 1//c , where W ±1 and A ±1 are amplitude parameters. 
These are the same solutions as in the textbook procedures. 

The point now is that the Riccati solution u p = i^- mentioned above is only the 
particular solution, i.e., u p l = — tan(cr/) and u p _ { =coth(c77) for K = ±1, respectively. 
The particular Riccati solutions are closely related to the common factorizations of the 
second-order linear differentail equations that are directly related to the well-known Darboux 
isospectral transformations [4]. Indeed, Eq. (6) can be written 

w" -c(-Kc)w = (7) 

and also in factorized form using Eq. (6) one gets (D^ = 4^-) 

(D n + cu p ) (Djj - cu p ) w = w - c(u p + cu 2 ) w = . (8) 

To fix the ideas, we shall call Eq. (8) the bosonic equation. On the other hand, the 
supersymmetric partner (or fermionic) equation of Eq. (8) will be 

(Djj - cu p ) (D,, + cu p ) w f = w f - c(-Up + cup)w f = w f - c • c K)SUSy w f = .(9) 
Thus, one can write 

, s ' 2 f c(l +2tan 2 cTl) if K = 1 

C,c,susy(77) = -Up+CU p 



c(-l +2coth 2 c7]) iffc = -l 

for the supersymmetric partner adiabatic index. The solutions w f are w f = cos ( c ^ +d -) an d 
w f = sinh ^^ for K = 1 and K = — 1 , respectively. 

Introducing the (quantum momentum) operator = — z'D^ we can write the fermionic 
equations as follows 

(-Pjj - icup) (P n - icup) w f = -P^w f - c(-iP r? u p + cu 2 ) w f , (10) 

whereas the bosonic case is 

(Pj, -icup) (-P n -icup) w b = -Pjw b -c(iP^u p + cUp)w b , (11) 
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(13) 



There is a more general factorization of the bosonic equation [5] 

(D^ +cu g ) (Dj, -cu g ) w g = Wg -c(u g + cug)wg = w g + kcc(t]; A)w g = ,(12) 
which is given in terms of the general Riccati solution u g (T]) 

u g (T];A) = u p (77) - [inCkft) + A)] = D„ [in L^L ^j 
and yields the one-parameter family of adiabatic indices c^rj; A) 

-kc k {t]-X) = c U 2(tj; A) + dUg ^ ;A) = -ice - -Dl [ln(I K (Tj) +A) 
5 dr/ c ' L 

4w r (T 7 )w / y (T]) 2w*(tj) 

C(I K (7])+A) + c(I (f (7 7 )+A) 2 ' V ; 

where Ik(tj) = /J 7 w^.(y)dy, if we think of a half line problem for which A is a positive 
integration constant thereby considered as a free parameter of the method.. 

All c K (rj ; A ) have the same supersymmetric partner index c KjSUS y (rj ) obtained by deleting 
the zero mode solution w K . They may be considered as intermediates between the initial 
constant index ?cc and the supersymmetric partner index c Ksusy (T]). From Eq. (12) one can 
infer the new parametric 'zero mode' solutions of the universe for the family of barotropic 
indices c K (r}\ A) as follows 

w 8 (.;A) = 7 #^^a 8 (,,A)=f^# T V. (15, 



Before closing this section, we recall an interesting point. Since what we have done here 
is to use the Darboux transformations at the level of cosmological evolutionary equations (i.e., 
equations of motion of the scale factor of the FRW cosmologies) a natural question is what 
is the effect of such transformations at the level of any Lagrangian of the cosmological fluid 
mechanics. The answer to this question has been already provided in the literature. Neto 
and Filho [6] have shown that in general the application of the Darboux transformations is 
equivalent to the addition of a total time derivative of a purely imaginary function to the 
Lagrangian and later, Samsonov [ ] using the coherent state approach confirmed their result. 

3 - Dirac-like formalism 

The Dirac equation in the susy nonrelativistic formalism has been discussed by Cooper et al 
[8] already in 1988. They showed that the Dirac equation with a Lorentz scalar potential is 
associated with a susy pair of Schroedinger Hamiltonians. This result has been used later by 
many other authors. Here we make an application to barotropic FRW cosmologies that we 
find not to be a trivial exercise except for the uncoupled 'zero-mass' case (subsection 3.1). 

3.1- Let's introduce now the following two Pauli matrices a = —\o y = — i ^ ? ^ J 

and j8 = a x = ( ^ q ^ and write a cosmological Dirac equation 

h frw w = + ax ( icUp )] W = o , (16) 

where W = ( ^ is a two component 'zero-mass' spinor. This is equivalent to the 

following decoupled equations 

— P^Wj +icu p Wj = (17) 
H-P^Wj + icu p w 2 = . (18) 



Solving these equations one gets Wj oc 1/cos(ct]) and w 2 oc cos(cr/) for K = 1 cosmologies 
and Wj oc l/sinh(c7]) and w 2 oc sinh(cr/) for K" = — 1 cosmologies. Thus, we obtain 



W 



Wo 



w, 



This shows that the matrix 'zero-mass' Dirac equation is equivalent to the two linear second- 
order differential equations for the bosonic and fermionic components. 

3.2- Consider now a "massive" Dirac equation 

h frw w = ^^p^ + CTx ( icUp + K )] W = KW ^ ( 19) 

where K may be considered the mass parameter of the Dirac spinor. Eq. (19) is equivalent to 
the following system of coupled equations 

— Ptj w j + (icup + K) w j = Kw 2 (20) 
Ptj w 2 + (icu p + K)w 2 = Kwj . (21) 

These two coupled first-order equations are equivalent to the following second order equations 
for the two spinor components 



-P>i-c 



i(^P r? -2K)u p + cu 



W; =0 



where the subindex / = 1,2. 



The fermionic spinor component can be found directly as solutions of 



D 2 w+ 



c 2 ( 1 + 2 tan 2 C7] ) + 2icKtan ct] 



w, 







and 



c 2 (-l +2coth 2 c7]) -2icKcothc7] 
whereas the bosonic components are solutions of 



w, 







for K = 1 



forK 



and 



D 2 w+ + 



D> 2 + 



c — 2icKtanc77 



w. 







-c +2icKcothc77 



w n 







for K — 1 



for K 



1 . 



(22) 

(23) 
(24) 

(25) 
(26) 



The solutions of the bosonic equations are expressed in terms of the Gauss hypergeometric 
functions 2 F l in the variables y = e lc7J and y = e c?7 , respectively 



1 

wJ = Ay p 2 Fi -2(P + i£ l) ; -o(P- i£ l)' 1 -P ; -y 



1 

2' 



+ 



B y P 2 F i (p - iq) , ^ (p + iq) , 1 + p; ~t 



(27) 



and 



w 



c(-i)-*y-Vi 



- ^(r + is), -^(r — is), 1 -ir;y 2 



+ 



D(-l)iy r 2 F 1 [i(r-s),^(r + s),l+ir;y 2 



(28) 
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respectively. The parameters are the following: p=(— 1 — ^)2, q= (1 — 

r = ( — 1 — i 2 ^) 2 , s = ( — 1 +i^f)K whereas A, B, C, D are superposition constants. 

It is not necessary to try to find the general fermionic solutions through the analysis 
of their differential equations (23) and (24) because they are related in a known way to the 
bosonic solutions. [9] The general fermionic solutions can be obtained easily if one argues 
that the particular fermionic zero mode is the inverse of a particular bosonic zero mode and 
constructing the other independent zero mode solution as in textbooks. Thus 



± _l+kf[w±] 2 dz 



w 



w 



± 



(29) 



where k is an arbitrary constant. 



3.3- The most general case in this scheme is to consider the following matrix Dirac-like 
equation 



-i 

1 



1 

1 



icUp + Kj 
icUp + K, 



Wj 

w n 



Kj 
K 



w- 



Proceeding as in 3.2 one finds the coupled system of first-order differential equations 



+ icug + K 2 



w- 



P,j + icUp + Kj 



w, 



KjWj 



K 9 w 2 



and the equivalent second-order differential equations 



-V w <+ 



ic(u p -Ug) + (K 1 -K 2 ) 



(30) 

(31) 
(32) 



^(iPjjUj. + KjUg + KjUp) -c 2 u p Ug 



w, = 



(33) 



where the subindex i = 1,2, and Uj and u 2 correspond to u p and u g , respectively. In the 
notation this equation reads 



cAupg-iAK) DjjW.-I- 



+ 



c(±D,jU f + (iKjUg + KjUp)) -c 2 u p u g 



w. = . 



Under the gauge transformation 

1 n 



one gets 



where 



w . = z . eX pi -- 



cAupg-iAK) dx \ =z i (ri) 



g jirjAK 



(1 K + X)2 

Pn 2 z i + fii(T?)z i = 0, or 0^ + 2,(77)^ = 0, 
Q f (n) = c(±D )? u. + (iK 1 Ug + K 2 Up))-c 2 UpU ? 



(34) 



(35) 



(36) 
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cAupg — iAK) 



(37) 



for / = 1 , 2, respectively. Q t are complicated 'potential' functions and we were not able to find 
analytical solutions of Eq. (36). 

The corresponding Dirac spinor is of the following form 



where w g (A,K 2 ) is given by Eq. (15) for Kj = K 2 = 0. For A — > °° one obtains W(A,0,0) — > 
W. In addition, u g — > u p and for Kj = K 2 = K one gets the particular case in 3.2. 

4 - Conclusions 

We come now to the interpretation of the mathematical results that we displayed in the 
previous sections. An examination of the formulas (23-26) and (37) show that the parameters 
K introduce an imaginary part in the adiabatic index of the cosmological fluid. Thus, the 
supersymmetric techniques presented in this research letter are a particular way to consider 
dissipation and instabilities in the ideal case of barotropic FRW cosmologies. More general 
scale factors of barotropic FRW universes incorporating a well-defined type of dissipation can 
be obtained from the 'zero-modes' Wj 2 ± by means of the relation a ~ w 1 / c ( T ' ;K ). The indices 
c(rj ; K) are redefined adiabatic indices that can be infered from the formulas (23-26) and (37), 
respectively. 
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